%per bus kijken wat de spanning doet bij een uitval van een lijn.
%nu nog zonder vermogenvariatie

clear all
tic
IEEE_14_bus_original;
einde_interval = size(Line.con,1);
aantal_lijnen = einde_interval;
teller = 0;
tijd = zeros(200,einde_interval);

alle_spanningen = zeros(200*einde_interval,14);

for k=1:einde_interval;
    
    filename = sprintf('lijn_uit%03d',k);
    initpsat;
    Settings.plottype = 2;
    Settings.freq = 60;
    Settings.absvalues= 'on';
    runpsat(filename,'pf');
    Settings.fixt = 1;
    Settings.tstep = 1;%0.05; %TRAAAAAAAAAAAAAG
    
    Settings.freq = 60;
    Settings.absvalues= 'on';
    Settings.plottype = 2;
    runpsat('td');
    
    %begint vanaf index 0... (STOM!!!)
   v = genvarname('voltages_fout_bus', who);
   eval([v ' = Varout.vars']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
   lengte = size(Varout.vars,1);
   
   %spanning per tijdsstap per bus
   alle_spanningen(1+200*teller:200*teller+lengte,:) = Varout.vars(:,1:14);
   
   t = genvarname('tijd_fout_bus', who);
   eval([t ' = Varout.t']);%(:,(size(Bus.con,1)+1):2*size(Bus.con,1))']);
   
   %tijd per stap per lijn
   tijd(1:lengte,k) = Varout.t;
   closepsat;
   teller = teller+1;
   
   save lijn_uit tijd alle_spanningen aantal_lijnen einde_interval;
end


% 
% 
% 
% 
% 
% 
% 
% 
% %waar bevinden zich de dips?
% dips_bussen_belasting = zeros(100,6);
% teller = 0;
% 
% for l = 1:aantal_lijnen%((aantal_lijnen)*7)
%     for ll = 1:aantal_bussen
%         if voltage_k_percent(l,ll)<=0.9 && voltage_k_percent(l,ll)>=0.01
%             teller = teller+1;
%             %1e index van dips geeft geval weer, 2e index de bus
%             dips_bussen_belasting(teller,1) = l;
%             dips_bussen_belasting(teller,2) = ll;
%             dips_bussen_belasting(teller,4) = voltage_k_percent(l,ll);
%             
%         end
%     end
% end
% aantal_dips = teller;
% %nullen verwijderen
% for l=1:size(dips,1)
%     if dips(l,1) == 0
%         
%         dips(l:size(dips,1), :) = [];
%         break;
%     end
% end
% 
% %er zijn grenzen opgegeven van min en max voltage per belasting (0.9-1.1)
% %checken of er een bus te laag of te hoog gaat:
% lijst_getroffen_bussen = zeros(100,2);
% teller = 1;
% for l=1:size(dips)
%     for k = 1:size(PQ.con)
%     if dips(l,2)==PQ.con(k,1)
%         lijst_getroffen_bussen(teller,2) = dips(l,2);
%         lijst_getroffen_bussen(teller,1) = dips(l,1);
%         teller = teller + 1;
%     end
%     end
% end
% 
% %overzicht van meest gevoelige bussen
% gevoeligheid_bussen = zeros(30,1);
% for l = 1:size(dips,1)
%     bus = dips(l,2);
%     gevoeligheid_bussen(bus,1) = gevoeligheid_bussen(bus,1) + 1;
%     
% end
% 
% toc;
% 
% 
% save lijn_uit dips_bussen_belasting aantal_dips 
% 
% 

